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FINITE  TEMPERATURE  STABIUZATION  OF  THE  GRADIENT 
DRIFT  INSTABILITY  IN  BARIUM  CLOUDS 

I.  INTRODUCTION 

'  The  study  of  the  dynamic  evolution  of  artificial  plasma  clouds  (e.g., 
barium  clouds)  in  the  earth's  Ionosphere  and  magnetosphere  has  been  an 
active  area  of  research  for  more  than  two  decades.  The  Initial  motivation 


for  these  active  space  experiments  was  to  use  the  plasma  cloud  as  a 
diagnostic  to  determine  ambient  plasma  conditions  (e.g.,  electric  fields, 
neutral  winds).  was  se^  discovered  that  the  Interaction  of  an 
artificial  plasma  cloud  and  the  near-earth  space  plasma  was  very  complex, 
and  that  plasma  clouds  undergo  a  complicated  and  dynamic  evolution.  One  of 
the  prominent  features  of  cloud  evolution  In  the  Ionosphere  Is  the 
development  of  field-aligned  strlatlons  on  the  cloud's  steepened 
"backside”.  These  strlatlons  have  been  attributed  to  the  onset  of  the 
gradient  drift  Instability^  (Llnsou  ~and~vr6fkfflaTt; — This  instability 
can  occur  In  a  weakly  colllslonal  plasma  which  contains  a  density  gradient 
and  a  neutral  wind  (or  ambient  electric  field),  (^mon; — t‘9b3r'?oh,  1961) .  r. 


51monr-t^b^5oh,"7r967)>  .  r  _ 
comniitatlonal  research  has  w 


A  substantial  amount  of  theoretical  and  computational  research  has  ^ 

been  Invested  In  understanding  the  linear  and  nonlinear  development  of  the“ 
gradient  drift  Instability  and  its  relevance  to  plasma  cloud  structure 
(Volk  and  Haerendel,  1971;  Perkins  et  al.,  1973;  Zabusky  et  al.,  1973;  j 

Shlau  and  Simon,  1972;  Perkins  and  Doles,  1975;  Scannapleco  et  al.,  1976; 

Chaturvedl  and  Ossakow,  1979;  Kesklnen  et  al.,  1980;  McDonald  et  al.,  1980; 

McDonald  et  al.,  1981;  Hubs  et  al.,  1983).  The  bulk  of  analyses  to  date 
have  neglected  the  effects  of  plasma  dynamics  parallel  to  the  ambient 
magnetic  field  Bq  (l.e.,  considered  perturbations  only  In  the  plane 

1 

transverse  to  ^)  or  have  Incorporated  parallel  effects  In  a  crude  way. 
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Francis  and  Perkins  (1975),  for  example,  assume  that  the  ambient  and 
perturbed  potentials  map  uniformly  along  Bq  thereby  connecting  the  cloud 
dynamics  to  that  of  the  conducting  background  at  different  altitudes.  In 
this  model  finite  thermal  effects  have  a  stabilizing  Influence  on  short 
wavelength  modes.  However,  several  studies  have  been  performed  which 
attempt  to  Include  parallel  dynamics  self-conslstently  In  the  stability 
analysis  (Goldman  et  al.,  1976;  Sperling  and  Glassman,  1984;  Sperling  et 
al.,  1984).  These  studies  Incorporate  the  parallel  length  of  the  cloud 
along  the  ambient  field  Into  the  stability  analysis  and  have  shown  that 
parallel  effects  can  play  an  Important  role  In  the  development  of  the 
gradient  drift  Instability  as  It  relates  to  Ionospheric  plasma  clouds. 

In  particular,  Sperling  et  al.  (1984)  have  recently  shown  that  when 
the  finite  size  of  the  cloud  along  Bq  Is  self-conslstently  incorporated  in 
the  linear  stability  analysis,  the  long  wavelength  modes  tend  to  have  much 
smaller  growth  rates  than  short  wavelength  modes.  This  reduction  In  growth 
due  to  finite  cloud  length  occurs  because  the  Integrated  Pedersen 
conductivity  of  the  background  plasma  over  the  extent  of  the  mode  along  the 
field  becomes  greater  than  the  cloud.  The  distance  the  mode  extends  along 
the  field  Is  proportional  to  the  perpendicular  wavelength  of  the  mode  so 
that  long  wavelength  modes  are  more  strongly  affected.  However,  this 
analysis  neglected  the  effects  of  finite  temperature.  Sperling  and 
Glassman  (1984)  Included  finite  temperature  in  the  analysis  and  found  that 
short  wavelength  modes  became  propagating  rather  than  purely  growing 
modes.  Moreover,  for  sufficiently  short  wavelengths  there  was  some 
evidence  that  the  modes  were  be  completely  stabilized.  A  drawback  of  this 
work  though  is  that  a  relatively  complex,  second-order  differential 
equation  Is  solved  numerically  to  obtain  results.  The  underlying  physics 
of  the  stabilization  mechanism  is  therefore  somewhat  obscure  and  not 


addressed  In  the  paper.  Nonetheless,  It  Is  evident  that  finite  temperature 
and  finite  cloud  length  effects  can  impact  the  development  of  the  gradient 

drift  instability. 

The  purpose  of  this  paper  is  to  present  a  linear  stability  analysis  of 
the  gradient  drift  instability  which  Incorporates  both  finite  temperature 
and  finite  parallel  length  effects.  A  simple  plasma  model  is  used  (similar 
to  the  one  used  in  Sperling  et  al.  (1984))  which  permits  an  analytical 
solution  to  the  dispersion  equation.  It  is  found  that  the  short  wavelength 
modes  are  stabilized  by  redistribution  of  electrons  parallel  to  Bq  (i.e., 
parallel  electron  diffusion  or  parallel  electron  response  to  the  perturbed 
fields)  which  neutralize  the  charge  Imbalance  set  up  by  the  Instability. 
An  analytical  expression  for  the  critical  wavenumber  for  the  stabilization 
of  instability  is  given,  as  well  as  numerical  results. 

The  organization  of  the  paper  is  as  follows.  In  the  next  section  we 
derive  a  set  of  general  nonlinear  equations  which  describe  the  evolution  of 
a  three  dimensional  plasma  cloud.  In  Section  III  we  discuss  the 
equilibrium  to  be  used  in  the  instability  analysis.  In  Section  IV  we 
derive  a  general  dispersion  equation,  and  present  both  analytical  and 
numerical  results.  Finally,  in  the  last  section  we  summarize  our  findings 
and  discuss  the  relevance  of  our  results  to  cloud  structure. 

II.  GENERAL  EQUATIONS 

We  first  derive  a  set  of  nonlinear  equations  to  describe  the  evolution 
of  a  warm  plasma  cloud  in  a  uniform  magnetic  field  B  =  Bq  z  with  a 

A 

background  neutral  wind  *  ts®®  Fig.  la].  For  simplicity  we 

consider  only  low  frequency  3/3t  «.  v ^  motion  of  the  cloud  and  take  the 
electron  collisions  to  be  sufficiently  weak  so  that  1  but 

allow  to  be  arbitrary.  The  collision  frequency  and  gyrof requency  of 


the  a  species  are  given  by  v  and  (2  ,  respectively.  Both  electrons  and 

n  Gi  I 

Ions  are  taken  to  be  warm  and  for  simplicity  we  consider  the  Isothermal 
limit.  In  this  case  the  fundamental  equations  of  our  analysis  are 
continuity,  momentum  transfer,  charge  neutrality  and  Ampere's  law; 
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where  and  are  the  fluid  velocity  and  temperature  of  species  a, 
respectively,  Is  the  a  specie-neutral  collision  frequency,  Is  the 

electron-ion  collision  frequency,  Is  the  Ion-electron  collision 

frequency,  and  a  refers  to  electrons  (e)  or  Ions  (l).  We  take  the  electric 
and  magnetic  fields  to  be  represented  by  potentials  as 

,  3A  - 

E  -  -  7^  -  -  -r;^  2  (6) 

~  c  3t 


and 


B  ••  B.S  z  +  7A  X  z 


(7) 


where  is  the  electrostatic  potential  and  A_  Is  the  vector  potential 

z 

associated  with  the  magnetic  field  produced  by  the  self-consistent  plasma 
currents.  We  consider  only  since  assume  1 7A^  x  z|  «  Bq. 

The  electron  cross-field  motion  Is  given  by 


V  7,ij»  X  z  +  ——  7,n  X  z 

~el  B  1^  eBn  1 


while  the  parallel  motion  is  given  by 


V,  =  -  [eE ,  +  T  V  ln(n)  +  f  (v  /v  )T7  tn(n)]/fm  (v  +  fv  .)1 
le  ‘■lei  ^  el  en-*  I  ■'  e  e  ei^- 


1  3A 


A  A 

where  v  =  v  +  v  ,  7  »  b  •  7,  b  -  B/|B^|,  and  T  =*  T  +  T . . 

6  6X  6I&  H  6  X 

The  assumption  that  1  while  ~  1  requires  that 

f  =  m  V  /n.v .  «  1 . 

e  en  1  In 

The  ion  cross-field  motion  is  given  by 
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2  2  — 

where  5.  *  fl  +  ,  and  the  parallel  motion  Is  given  by 

1  in  1  ^ 


''li  ’  HeE,  +  TJ,  in(n)]v^^  -  v^T7,  in(n) }/ [m^v^^[v^  +  fv^^)], 


In  (11)  we  have  included  both  the  Pedersen  and  Hall  responses  to  the 
electric  field,  neutral  wind  and  pressure. 

Substituting  (8)-(12)  into  (1),  (4)  and  (5)  we  find  that 


4^-:§7(|.xz  •Vn  +  V  T^A  -  D.V.n)  »  0, 

dt  B  1  '4 ire  1  z  II  I 


(13) 


2 
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7n 
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dA 


n2,  4ir  ^  1  z^ 

l^z  *  CP  ^^ll'**  c  dt 


(15) 


2 

where  *  m^u^/ne  is  the  parallel  resistivity,  0^  “ 

Di  *  5, fv,  /n. IcT/eB  are  the  parallel  and  perpendicular  ion  transport 
1  i ^  in  1  ^ 

coefficients,  T  *  +  T^,  and 
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(16) 


dt 
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7. 


(18) 


Equation  (13)  is  the  electron  continuity  equation,  (14)  is  the  charge 
conservation  equation  and  (15)  is  Ampere's  law.  Note  that  the  electron 
pressure  has  been  absorbed  into  (>  in  (16)  and  that  terras  of  order  f  «  1 
have  been  discarded  compared  with  those  of  order  unity.  Equations  (13)- 
(18)  constitute  a  complete  description  of  the  evolution  of  a  three- 
dimensional  warm  plasma  cloud. 
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III.  EQUILIBRIUM 

We  will  consider  the  linear  stability  of  a  two-dimensional  cloud  which 

is  localized  both  along  and  across  the  magnetic  field  Bq:  n^,  =  nj.(x,  z)  with 

n^,  ^  0  for  |x|  <  x^  and  |zl  <  z^.  The  background  plasma  is  taken  to  be 

uniform  throughout  the  region  I z I  <  z,  between  two  insulating  plates  at  z  = 

D 

+  z,  .  The  location  of  the  plates,  enables  us  to  control  the  ratio  of  the 
total  magnetlc-fleld-line  Integrated  Pedersen  conductivity  of  the  cloud  to 
that  of  the  background,  ^cZ^^/Cn^z^  +  n^z^j).  This  ratio  is  an  important 
parameter  of  the  equilibrium  configuration. 

The  equations  describing  the  two  dimensional  plasma  cloud  are  given  by 


il-D  i-  -  0 

3t  ”l  3^2  4..  3.  3,2 


2 

,  c  ^in  3  3  a  j.  n  4.  ^  ^ 

^i  B  3x  3x  ®1  g^2-  4  ire  3z  g^2 


0, 


3^A 


3A 


3.2  cn. 
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(19) 

(20) 

(21) 


When  there  are  no  density  gradients  in  the  x  direction,  the  solution  to 
these  equations  is  (|)  =  “  (B/c)  V^x  so  that 


(22) 

with  A^  =  0,  l.e.,  the  potential  adjusts  itself  so  that  the  electrons  are 

in  force  balance  along  the  magnetic  field  as  the  cloud  diffuses  along 
In  writing  the  solution  for  in  (22)  we  have  required  <p  0  for  |z|  ^ 
since  we  have  assumed  that  only  the  charge  accumulations  associated  with 
the  motion  of  the  plasma  cloud  itself  are  responsible  for  the  development 
of  (>.  There  is  no  charge  at  z  =  +  z^.  The  density  satisfies 
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0 


(23) 


3n/3t  -  D.3^n/3z^ 
If 


50  that  n  diffuses  along  In  the  opposite  limit,  where  there  are  no 

iensity  gradients  in  the  z  direction  ( 3/ 3z  =  0),  A^,  =  0  and 
(t  *  -  (T/e)!ln(n)  -  (B/c)(v.  /fJ.  )V  dx(n,/n)  so  that 


♦  -  -  ^  ^[”^1  *  1  >  v„  /“d,f.  -  ::i). 


"b  ■  = 


(24) 


where  we  have  required  0  as  | x j  In  Eq.  (24)  the  potential  adjusts 

itself  so  that  the  perpendicular  ion  pressure  gradient  is  balanced  by  the 
electrostatic  field.  In  comparing  the  expressions  for  i)>  in  (22)  and  (24), 
it  is  important  to  note  that  the  electron  and  ion  pressures  push  the 
potential  in  the  opposite  direction  (compare  the  signs  of  the  terms 
proportional  to  Tg  and  T^). 

We  now  return  to  the  more  general  two  dimensional  equations  in  (19)- 

(21).  These  equations  describe  three  .  time  scales:  the  resistive  flux 

2  2  2 
diffusion  time,  T  =  4ir  x  /p  c  ;  the  parallel  diffusion  time,  x,  =  z  /D„; 

r  c  e  II  c  If 

2 

and  the  perpendicular  diffusion  time,  We  assume  that  the  flux 

diffusion  time  is  the  shortest  time  scale  so  that  inductive  effects  are  not 
important  in  the  equilibrium,  l.e.,  in  Eq.  (21)  3A^/3t  «  c  3(i)/3z.  The 
equations  then  simplify  to 


2  2 

3n  _3n  _3n  ,c3  3(t  ^ 
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(25) 
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[.  CONCLUDING  REMARKS 


We  have  presented  a  relatively  simple  analysis  of  the  gradient  drift 
nstability  in  barium  clouds  which  includes  both  finite  temperature  and 
inite  parallel  length  effects.  We  have  derived  a  general  set  of  nonlinear 
quations  which  describe  the  evolution  of  a  three  dimensional  plasma  cloud 
n  the  ionosphere.  We  then  investigate  the  stability  of  a  two  dimensional 
lasma  cloud  in  which  the  density  varies  in  the  direction  of  the  neutral 
rind  and  along  the  ambient  magnetic  field  Bq.  By  modeling  the  density 
'arlation  along  Bq  as  a  waterbag,  we  are  able  to  obtain  an  analytic 
llsperslon  relation  for  the  gradient  drift  instability.  There  is  no  static 
squlllbrlum  since  the  cloud  diffuses  both  perpendicular  and  parallel 
:o  Bq  so  that  we  restrict  our  analysis  to  time  scales  short  compared  to  the 
:loud  diffusion  time  scales. 

In  recent  work  Sperling  et  al.  (1984)  found  that  the  finite  length  of 
the  cloud  suppressed  the  growth  of  long  wavelength  modes  because  of  "good" 
coupling  to  the  background  plasma.  However,  the  growth  rate  of  short 
wavelength  modes  is  unaffected  by  the  finite  size  of  the  cloud  since  they 
do  not  couple  to  the  background  in  a  zero  temperature  plasma.  In  contrast 
to  this  result,  we  find  that  in  a  finite  temperature  plasma  the  short 
wavelength  modes  do  couple  to  the  background  plasma  and  that  for 
sufficiently  short  wavelengths  the  modes  are  completely  stabilized. 
Stabilization  results  as  the  electrons  redistribute  parallel  to  B^  and 
neutralize  the  charge  imbalance  set  up  by  the  instability.  The  parallel 
electron  motion  results  from  diffusion  and/or  the  response  to  the  perturbed 
fields.  We  have  also  derived  a  simple  analytic  expression  (see  (54))  for 
the  wavenumber  corresponding  to  marginal  stability  (i.e.,  y  ®  0)*  Our 
results  are  consistent  with  previous  numerical  computations  of  Sperling  and 
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V. 


NUMERICAL  RESULTS 


We  now  present  quantitative  results  for  the  wave  frequency  of  the 

gradient  drift  instability  by  solving  (43)  numerically.  In  Fig.  2  we  plot 

y/y«  (solid  curve)  and  <»>  /Yf>  (dashed  curve)  vs.  k  p  where  aj  =  cu  +  iy  is 

the  wave  frequency.  We  consider  the  following  typical  parameters  for  a 

barium  cloud  at  180  km:  n  /n.  ■  10.0,  T  “T  *0.1  eV,  m^  *  16  m  (0^ 

CD  el  H 

background),  *  10  km,  B  *  0.3  G,  a  “  10^,  *  100  m^/sec,  =  '^eff  * 

20  m/sec  and  L„  =  1  km.  From  these  values  we  note  that  c  ®  10^  m/sec  and 

ll  S 

p  =  6  m.  The  main  features  of  Fig.  2  are  described  as  follows.  First, 
s 

the  growth  rate  has  a  maximum  value  at  k  p  “  k  _^p  *  0.08.  For  k  <  k_ 

°  y  s  ym  s  y  y™ 

the  growth  rate  decreases  because  of  coupling  to  the  background  plasma  as 

described  in  Sperling  et  al.  (1984).  For  ky  >  ky^  the  growth  rate 

decreases  rapidly  and  becomes  stable  (i.e.,  y  <  0).  This  is  due  to 

parallel  electron  motion  as  described  in  the  preceding  section  and  is  the 

dominant  finite  temperature  effect.  Second,  the  critical  wavenumber  for 

stabilization  as  given  by  (54)  is  denoted  by  the  arrow  along 

the  k  p  axis  fk  p  a  0.26).  It  is  seen  that  (54)  gives  a  very  good 

y  s  y  3 

approximation  to  the  critical  wavenumber  obtained  numerically 

(kyp^  “  0.28),  For  the  parameters  used  the  critical  wavelength  is  given 

by  X  a  135  m.  Third,  the  real  frequency  is  linear  in  k„  and  is 

c  •  J 

proportional  to  the  diamagnetic  drift  velocity  Vj  =  (cT/eB) fr//n^) . 

The  results  shown  in  Fig.  2  appear  to  be  in  qualitative  agreement  with 
the  numerical  work  of  Sperling  and  Classman  (1984).  They  consider  a  plasma 
cloud  with  a  Gaussian  distribution  along  the  magnetic  field  rather  than  a 


waterbag  model.  The  shape  of  the  curve  of  y  vs  k^  shown  in  Fig.  2  is 

similar  to  corresponding  curves  presented  in  their  report.  Moreover,  they 

also  find  that  lu  ®  k  for  large  k„. 

r  y  y 
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with 


n 

% 


iVcri 

oD  , 


(55) 


where  L„  =  -  fn  "/n  1  ^  and  p  *  c  /fi,  is  the  Larmor  radius  based  on  the 
n  cc  s  si 

sound  speed  *  (T/m.  Note  from  (54)  that  in  the  limit  T  0,  then 

kj.  "  so  the  mode  is  not  stabilized  for  a  zero  temperature  plasma 
[Sperling  et  al.  (1984)].  Equation  (54)  also  appears  to  imply  that  k^ 

becomes  very  large  as  the  integrated  cloud  density  n^  Increases.  This 
result  is  misleading.  In  a  cloud  which  is  finite  in  the  y  direction,  the 

neutral  wind  velocity  in  (54)  should  be  replaced  by  the  effective  slip 

velocity  of  the  cloud  and  the  neutral  wind.  For  a  circular  cloud,  the 

effective  velocity  is  given  by  [Ossakow  and  Chaturvedi,  1978] 


eff 


2V  n.  2, 
n  b  b 


/(^Vb  Vc^* 


(56) 


In  the  limit  of  very  large  n^^c*  yields  the  critical  wavelength 


2  2 
k  L 
c  n 


2 

4V  z. 
n  b 

2D^ 


(1 


/  2„2 

4z,  V 

+  ->■ 

2  2-’’ 
ap  c 
s  s 


(57) 


which  is  independent  of  nj.Zj,.  Thus,  short  wavelength  modes  are  stable  even 
in  large  ionospheric  clouds. 
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Two  significant  features  of  (52)  are  first  that  Yq  does  not  enter  the 
equation,  and  second  that  there  appears  to  be  an  instability  even  in  the 
absence  of  Yq*  However,  the  unstable  root  is  invalid.  In  order  for  to 
remain  bounded  as  |z|  becomes  large,  we  require  Re(k^j)  >  0.  Thus, 
for  Y  “  im*  we  find 


(ui*/aD^) 


1/2  iir/4 


The  root  with  the  opposite  sign  must  be  discarded.  Inserting  this 
expression  for  k^j  into  (46),  we  obtain  the  damped  root  in  (52).  Thus,  when 
thermal  effects  are  retained,  the  gradient-drift  mode  is  stable  for  large 


The  stabilization  mechanism  can  be  understood  in  terms  of  our  earlier 
discussion  of  the  local  dispersion  relation.  When  is  discarded  in  (45) 
and  (46),  the  parallel  wavevector  k^  remains  real  (the  mode  is  evanescent 
in  the  outer  region)  and  the  instability  survives  even  as  ky  •.  However, 
when  is  retained,  the  parallel  wavevector  k^  becomes  complex  (the  mode 
is  evanescent  and  propagating  in  the  outer  region)  and  the  mode  is  stable 
for  large  ky.  In  other  words,  cu*  causes  the  mode  to  become  oscillatory 
along  ^  and  the  electrons  can  now  move  along  rather  than  across  B  to 
neutralize  the  ions;  the  mode  is  therefore  stable. 

To  obtain  an  expression  for  the  critical  wavelength  k^,  where  the  mode 
becomes  stable,  we  take  Y  “  “  iu^  and  solve  (46)  directly.  The  critical 
wavelength  is  given  by 


4  [(1 


,22 
2n  z 
c  c 


V  , 
n 

oD  ,L  ^ 
1  n 


,  2  2  2  2 
4n  z  pc 
c  c  3  3  1 

2  ,,2.  2 
n,  oD  ,L 
b  In 
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.  -^2  Vr  - 

[Yq  -  yJ  -  a  -j-J  Y, 


2  2  2 

and  Is  independent  of  the  plasma  temperature.  For  z  »  oD  n, /n 

c  r  D  c  u 


T  -  Yr 


2  2  2 

while  in  the  opposite  limit  z  <<  oD  n, /n  y^. 

^  c  r  b  c  0’ 

2  2  2 

- 

on^r 

These  results  have  been  obtained  previously  [Sperling  et  al.  (198A)]. 


2  _ 

When  ky  becomes  very  large,  *  y/oDj^  and  the  dispersion  relation 


becomes 


Cr  -  -  Yg 


Y. 

n  z 
c  c 


For  large  ky,  we  have  »  Yg  since  ui^  «  k^  and  therefore  to  lowest  order 


y  *  i"*. 


the  mode  simply  propagates  at  the  diamagnetic  frequency.  To  investigate 
whether  the  mode  is  stable  or  unstable  in  this  limit,  corrections  to  this 
result  must  be  calculated.  These  corrections  arise  from  the  right  side  of 
(50)  and  yield  the  eigenvalue 


r  ^  j.  °  r  iiT/4 

c  c 


18 


where 


A 

k 


c  c  Yq  -  Y  +  loj* 


In  Sperling  et  al.  (1984),  the  dispersion  relation  was  carefully  solved  In 
all  possible  limits  of  the  various  parameters.  In  order  to  clearly 
emphasize  the  role  of  finite  temperature  in  stabilizing  short  wavelength 
modes,  we  will  focus  our  analysis  primarily  on  modes  with  growth  rates  less 
than  Yq  for  which  kj^/k^,  «  1.  In  this  limit  (43)  becomes 


k  z  -  k, 

C  C  D 


Vc^  -  la.*  -  Yg)  +  n^^k;^  Y 


In  the  form  shown  in  (46),  the  nature  of  the  dispersion  relation  can  be 
simply  understood.  The  first  (second)  term  is  the  magnetic  field  line 
integrated  contribution  from  the  cloud  (background).  The  background 
contribution  arises  only  from  the  region  [zl  <  k^^  since  the  perturbed 
potential  is  small  outside  of  this  Interval.  The  result  In  (46)  can  be 

A' 

directly  obtained  by  integrating  (33)  for  along  z. 

The  dispersion  relation  in  (45)  or  (46)  can  be  solved  analytically  for 

2  _ 

various  values  of  ky.  In  the  limit  ky  0,  k^  *  y/oD^  and  the  dispersion 
relation  is  given  by 


17 


discontinuity  In  the  density  at  z  *  +  z_,  we  find  that  and  A  must  be 

A  ^2 

continuous.  For  the  even  ^  solution  *  O),  we  obtain  the  dispersion 

equation 


1c  z 
c  c 


,  k.  Y  +  k  D 


+  miT 


c  Y  +  k  D 

'  XT 


-  f-f  * 

y  rb 


(42) 


where  m  is  an  Integer.  The  dispersion  relation  for  the  odd  ^  mode  (41^  ■ 
0)  is  similar  to  (42)  except  tan“^  is  replaced  by  -  cot  The  result  In 
(42)  is  Identical  to  the  dispersion  relation  obtained  previously  by 

Sperling  et  al.  (1984)  except  that  the  expressions  for  k^j  and  k^.  given  In 
(40d)  and  (41d)  now  contain  thermal  effects  which  were  previously 

neglected. 

In  general,  the  dispersion  equation  (42)  has  an  infinite  number  of 

solutions  for  a  given  set  of  physical  parameters,  corresponding  to 

elgenmodes  with  an  increasing  number  of  nodes  (m)  along  z.  The  dispersion 

equation  (42)  can  be  solved  numerically  for  arbitrary  values  of  the 

background  and  cloud  parameters.  However,  to  gain  an  understanding  of  the 

general  scaling  of  the  growth  rate  Y  with  the  parallel  extent  of  the  cloud 

and  temperature.  It  is  useful  to  solve  (42)  analytically.  To  do  this  we 

make  a  number  of  simplifying  assumptions.  We  consider  only  the  lowest 

order  mode  (l.e.,  m  =  0  which  implies  0  <  k^z^  <  it/2;  it  Is  easily  shown 

that  this  mode  has  the  largest  growth  rate);  take  v.  /7l.  to  be  small  so 

in  1 

that  take  all  parameters  but  the  density  to  be  the  same  inside  and  outside 
of  the  cloud;  assume  n_  »  nv;  and  v  ,  »  v  so  that  the  resistive 

diffusion  coefficient  is  continuous  across  the  boundary  z  *  z^ 
(l.e.  ).  With  these  assumptions  the  dispersion  equation  becomes 


k  z 
c  c 


tan"^ (k^/k^) 


(43) 
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with 


(^b  Vrb)*bb  ■  V*b 


(40c) 


2  Y(Y+kjD^) 

c,  »  - ^ -  I 

°  a  D  (7  +  k'^D.)  b 
r  y  i-'  ' 


(40d) 


where  the  subscript  b  on  a  given  parameter  indicates  that  it  is  to  be 

evaluated  in  the  region  |zl  >  z^.  Note  that  we  have  assumed  k^jZ^j  »  1  and 

that  solutions  which  diverge  as  z  becomes  large  have  been  omitted  from 

1/2 

(40).  This  assumption  is  consistent  with  the  approximation  z^  «  o  x^ 
made  in  deriving  the  equilibrium  potential  as  long  as  k^jZ^j  »  1. 

In  the  region  |  z  |  <  z^  the  solutions  for  and  are 


A'  ■  A  ^  sin  (k  z )  -  A  ^  cos  (k  z ) 
z  zc  ^  c  •'  zc  ^  c  ■' 


(41a) 


<j>  -  cos  (k^z)  +  ((i^  sin  (k^z) 


(41b) 


(y  +  k^  D  )a^’^  =■  k  c*^ 
'  c  y  rc-^  zc  c  c 


(41c) 


aDr(7  +  k^Dj 


(41d) 


We  note  that  in  writing  (40d)  and  (41d)  we  have  assumed  y  »  k^V^  for 
simplicity. 

To  complete  the  dispersion  relation,  we  must  match  the  various  plane 
wave  solutions  at  z  ■  +  z^..  The  appropriate  matching  conditions  are 
obtained  from  (32)  and  (33).  Integrating  these  two  equations  across  the 
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.  .I*  .• 


where  Djj^  ■  5^^  )cTj^/eB.  In  the  limit  »  0,  n^  *  (ikyCn‘'/YB)^ 

and  the  usual  growth  rate  of  the  instability  is  obtained  by  equating  n  and 
n^.  In  this  case  the  density  perturbation  n^  and  the  electric  field 
perturbation,  -  Ik^^,  are  in  phase  and  the  usual  physical  description 

of  instability  applies.  However,  in  the  limit  of  large  k_  with  T_  *  0, 
n^  *  n^e^/T^,  l.e.,  the  electrons  are  adiabatic  and  neutralize  the  ions  by 
moving  along  B  rather  than  across  B.  In  this  case  n^  is  out  of  phase 
with  Ey  and  there  is  no  instability.  Vhen  »  0  and  k^  is  large,  the 
mechanism  is  somewhat  different  although  the  mode  is  also  stable.  In  this 

-yj  ^ 

case,  the  electrons  bunch  parallel  to  B,  and  n^  ■  Ck^/en^Y)^*  Por  k^  large 

2 

the  ions  can  only  neutralize  the  electrons  if  Y““kD,.  +ikVv.  /fJ. . 

y  11  y  n  in  1 

Namely,  the  electrons  bunch  parallel  to  B  and  the  ions  diffuse  across  B  to 
neutralize  the  charge  Imbalance.  Finally,  we  note  that  while  electron 
parallel  diffusion  or  ion  perpendicular  diffusion  stabilizes  the  mode  when 
the  perturbation  is  periodic  along  ^  when  the  spatial  dependence  of  the 
mode  is  exponential,  i.e.,  p(z)  ~  p  exp  (;^  ^2^^*  thermal  effects  are 

not  stabilizing.  The  distinction  between  a  periodic  and  a  nonperiodic 
solution  is  Important  in  Interpreting  the  dispersion  relation  which  is 
obtained  for  the  equilibrium  shown  in  Fig.  lb. 

We  now  solve  (32)  and  (33)  for  the  profile  given  by  (19).  The 

^  ^  II 

boundary  conditions  used  are  0  as  |z|  z^.  For  the  step  profile 

for  n(x,  z)  the  solutions  to  (32)  and  (33)  in  the  region  |z|  >  z^  can  be 
written  as  plane  waves 


\  *  ^zb 


(40a) 


%  -  K  exp  [-  k. |z|  ] 


(40b) 


''  e  e  k-D  +  k‘-V‘/v,  5, 
y  r  z  A  In  1 

The  Alfven  wave  has  a  stabilizing  Influence  on  the  instability  in  the 
electromagnetic  limit  but  thermal  effects  do  not  affect  the  growth  rate. 

Of  course,  the  gradient  drift  instability  is  not  stable  when  ^ 

Tq  as  might  be  construed  from  (35)  since  the  expression  for  y  given  in  (35) 

breaks  down  when  y  <  Thermal  effects  can  have  a  significant 

2 

Influence  on  the  growth  rate  when  y  <  (electrostatic  limit).  Namely, 


k  T/m  V  <  y., 
zee  0 


the  electrons  can  move  a  wavelength  along  B  during  the  growth  time  of  the 
instability,  and  the  mode  is  stable.  In  the  limit  of  k  *  in  (36), 


-  k^D 

y  J- 


so  that  the  mode  damps  at  the  ion  diffusion  rate. 

The  stabilization  mechanism  can  be  understood  by  examining  the 
electron  and  ion  density  perturbations  in  the  electrostatic  limit. 


(ik  cn'/B  +  k  /en  ) 
y  z  e 


y  +  k  T  /m  V 
z  e  e  e 


(39a) 


[ik  cn"/B  -  kjc(v  /n  )n/B] 


ik  V  V,  /n. 

y  n  in  1 


In  deriving  (32)  and  (33)  we  have  neglected  parallel  ion  diffusion  as 
discussed  previously  and  assumed  n^.  »  n^.  The  primes  denote  a  derivative 
with  respect  to  x.  The  Important  finite  temperature  effects  that  appear  in 
(32)  and  (33)  are  the  diamagnetic  drift  frequency  (u*),  perpendicular  ion 

A 

diffusion  and  modification  of  the  equilibrium  In  the  limit  Tj^ 

“  Tg  =•  0,  Eqs.  (32)  and  (33)  reduce  to  those  previously  derived  by  Sperling 

et  al.  (1984).  The  eigenvalue  Y  is  in  a  frame  of  reference  moving  with  the 

electron  fluid. 

Prior  to  solving  (32)  and  (33)  for  the  density  profile  shown  in  Fig. 
lb,  we  first  consider  a  cloud  of  infinite  extent  [z^  **■  *]  and  Fourier 
expand  in  the  direction  parallel  to  l.e.,  p(z)  p  exp  [ik^zj.  This 

allows  insight  into  the  Influence  of  finite  temperature  and  parallel 

d3mamics  on  the  instability.  The  local  dispersion  equation  is  given  by 

k^V^ 

(T  -  io**  -  Yo)(?  +  k^D^)  .  -  (y  -  IkyV^  +  k^Dj,  (34) 


I  /2  ^ 

where  *  B/(4imm^)  is  the  Alfven  velocity.  This  dispersion  relation 
Illustrates  the  coupling  between  the  gradient  drift  instability  and  the 
Alfven  wave  when  k^  ^0. 

For  simplicity,  we  neglect  the  terms  in  (34)  which  cause  the  mode  to 

propagate  k^V^^  +  0  and  Y  ♦  y)  and  note  that  »  Dj^  in  ionospheric 

applications.  The  growth  rate  can  be  easily  obtained  in  two  limits:  in 

2 

the  electromagnetic  limit  (y  »  k  Dj,)» 


A 

while  in  the  electrostatic  limit  (y  «  k^D  ) 

V  '  y  f' 


(35) 


When  the  electron  parallel  conductivity  Is  large,  k  Is  typically  quite 

fl 

small  so  that  these  Inequalities  are  easily  satisfied.  Finally,  the  time 
evolution  of  the  equilibrium  should  have  a  negligible  influence  on  the 
stability  analysis  providing 


IYT|||,  Iyt^I  «  1 

where  y  is  the  growth  (damping)  rate  of  the  mode  of  interest. 


IV.  LINEARIZED  EQUATIONS  AND  DISPERSION  EQUATION 

To  determine  the  influence  of  finite  temperature  and  parallel  dynamics 
on  the  gradient  drift  instability  we  linearize  (13)-(15)  using  the 
equilibrium  discussed  in  the  preceding  section.  We  assume  perturbed 
quantities  to  vary  as  p  ~  p(z)  exp  (yt  +  ik^y).  After  eliminating  n  alge¬ 
braically  from  (13)-(15),  we  obtain  two  coupled  differential  equations 

A  ^ 

for  and  A  : 

2 


Bz 


(32) 


(7  -  IW*  -  yo)i  -  -  oD^(7  -  ik  +  k^D^)  i 


3A 
_ 2 

3z 


(33) 


where  y  *  y  +  i(c/B)ky<frQ  ,  - 


(T  /e}  ta 
e  ' 

-1 

1  •  “*  *  - 
and 


V 

n 


V  + 
n 


n  V,  _  n  n  z 

(1  In  y  _  T  _c  j  c  c 

'^cfi.  n  en  nz  +  n.  z. 
i  c  c  c  b  b 
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2  2 

„3An  ,c3.3<(» 

'»  - o  ^1  — o - 0^  ^  An  =  0, 

®  3z^  3x^  i  B  9x  3x 


2  ^  2 

1  3A<^__^3An  ,c'^i3  .  3<b 


where  A(ti  «  (^.  The  equation  for  An  has  no  steady  state  solution  so  An  will 
generally  evolve  on  a  time  scale  of  order  Xj  ~  The  z  dependent 

"  _i  *■ 

potential  A(^  drives  an  equilibrium  current  3A(|»/3z  which  is  required 

so  that  7«J  =  0.  The  time  variation  of  the  equilibrium  and  the  equilibrium 
parallel  current  J^,  will  be  neglected  when  we  carry  out  the  linear 
stability  analysis. 

Specifically,  in  the  linear  stability  analysis  we  formally  take  the 
limit  D  -►  0  and  choose 

R 

n^(x)  +  n^  |z|  <  z^ 

n(x,z)  =»  (31) 

”b  %  >  1*1  ^ 

as  shown  in  Fig.  lb.  The  perpendicular  diffusion  coefficient  is  retained 
in  the  stability  analysis.  The  step  function  model  for  the  density  is 
valid  as  long  as  the  parallel  wavelengths  of  the  modes  of  interest  are  much 
longer  than  the  actual  parallel  equilibrium  density  scale  length,  i.e., 

2  2 

k,  z^  «  1. 

I  c 


The  neglect  of  the  equilibrium  parallel  current  and  parallel  diffusion  can 
similarly  be  justified  for 


k.J  /ne  «  u 
n  z 


10 


Equation  (25)  governs  the  time  evolution  of  n  while  (26)  determines  (p.  For 
2  2 

z,  «  X  Q. /v  V.  ,  the  second  term  in  (26)  is  much  larger  than  the 

D  c  e  1  e  in 

remaining  terms  unless  3<^/3z  ■  0.  Thus,  4)  is  nearly  constant  along 
An  equation  for  can  then  be  obtained  by  averaging  (25)  and  (26) 
along  The  resulting  equations  are 


where  the  average  p  of  a  function  p(z)  is  defined  as 

_ 

p  *  /  dz  p(z)/2z^.  (29) 

“^b 

Thus,  the  Integrated  density  n  is  a  constant  in  time.  From  (28),  we  obtain 
an  expression  for  the  potential. 


V„  j"'  dx(l 


(30) 


The  first  and  second  terms  in  (30)  compete  in  forcing  the  potential  to 
balance  the  parallel  electron  pressure  and  ion  perpendicular  pressure, 
respectively.  In  the  limit  in  which  n  is  constant  along  Bq,  n  =  n  and  (ji  ~ 

-  tn  so  that  i  balances  the  ion  pressure.  In  opposite  limit, 

where  n  ■  and  )  ^  '^T^/e }  £n  fn/n^J  so  balances  the  electron 

pressure. 

Finally,  subtracting  the  field  line  average  of  (25)  and  (26)  from 
Itself,  we  obtain  an  equation  for  An  =  n  -  n  and  A(p  = 


9 


Giassman  (1984)  who  also  noted  the  Importance  of  finite  temperature  on  the 
stability  of  the  gradient  drift  mode. 

The  extent  to  which  the  work  described  here  may  have  relevance  to  the 
barium  cloud  "freezing"  phenomenon  deserves  some  discussion.  Briefly,  it 
is  observed  that  barium  ion  clouds  released  at  ionospheric  altitudes  will, 
through  the  nonlinear  evolution  of  the  gradient  drift  instability,  break 
into  smaller  pieces.  This  process,  known  as  bifurcation,  continues  until  a 
certain  minimum  transverse  dimension  is  reached,  at  which  point  further 
bifurcation  ceases.  This  frozen  scale  size  is  observed  to  be  approximately 
400  m.  Simple  two  dimensional  models  of  barium  cloud  evolution  fall  to 
explain  frozen  scale  sizes  this  large  [McDonald  et  al,  1981].  One 
explanation  which  has  been  advanced  to  explain  freezing  is  that  of  end¬ 
shorting  (Francis  and  Perkins,  1975;  Zalesak  et  al.,  1984]  which  takes  into 
account  the  distribution  of  plasma  along  magnetic  field  lines,  finite 
plasma  temperature,  and  parallel  currents,  as  we  have  done  here.  However 
end-shorting  depends  on  the  "mapping"  of  transverse  electric  fields  over 
relatively  large  distances  along  ^  this  becomes  increasingly  more 
difficult  as  k  Increases.  Indeed,  as  shown  here,  and  in  Sperling  et  al. 
(1984),  end-shorting  fails  to  provide  stabilization  of  the  very  high  k 
modes  (at  least  for  the  simple  profiles  of  n^.,  n^,  and  v  considered).  If 
nothing  else,  the  present  work  shows  the  existence  of  a  viable 
stabilization  mechanism  (i.e.,  finite  temperature)  in  the  gap  missed  by  the 
end-shorting  mechanism.  It  should  also  be  pointed  out,  however,  that  we 
are  not  that  far  removed  (135  m  stabilization  wavelength)  from  offering  a 
mechanism  that  by  itself  may  stabilize  a  400  m  diameter  barium  striatlon. 
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Given  the  uncertainties  in  applying  a  linear  theory  to  explain  a  complex 
nonlinear  phenomena,  the  stabilization  of  the  gradient  drift  Instability  by 
parallel  electron  diffusion  must  be  considered  as  a  strong  candidate  to 
I  explain  striation  freezing. 

I 
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555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  O.B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
01CY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOGSTON 
01CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Uilshipe  BLvd. 

Santa  Monica,  Calif  90401 
OICY  ATTN  C.B.  GABBARD 
OICY  ATTN  R.  LELEVIER 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 

■GENERAL  ELECTRIC  COMPANY  . 

SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
OICY  ATTN  M.H.  BORTNER 


BERKELEY  RESEARCH  ASSOCIATES,  INC.  SPACE  SCI  LAB 

P.O.  BOX  983 

BERKELEY,  CA  94701  GENERAL  ELECTRIC  COMPANY 

OICY  ATTN  J.  WORKMAN  P.O.  BOX  1122 

OICY  ATTN  C.  PRETTIE  SYRACUSE,  NY  13201 

OICY  ATTN  S.  BRECHT  OICY  ATTN  F.  REIBERT 


BOEING  COMPANY, 
P.O.  BOX  3707 

THE 

SEATTLE, 

WA 

98124 

OICY 

ATTN 

6. 

KEISTER 

01  CY 

ATTN 

D. 

MURRAY 

OICY 

ATTN 

G. 

HALL 

01  CY 

ATTN 

J. 

KENNEY 

GENERAL  ELECTRIC  TECH  SERVICES 
CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

OICY  ATTN  G.  HILLMAN 
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EOPHYSICAL  INSTITUTE 
NIVERSITY  OF  ALASKA 
AIRBANKS,  AK  99701 
(ALL  CLASS  ATTN;  SECURITY  OFFICER) 
OICY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
OICY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

;T£  SYLVANIA,  INC. 

LECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
7  A  STREET 
lEEDHAM,  MA  02194 
OICY  ATTN  DICK  STEINHOF 

(SS,  INC. 

!  ALFRED  CIRCLE 
JEDFORD,  HA  01730 

OICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 
01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
01CY  ATTN  J.M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFARD 
OICY  ATTN  JOEL  BENGSTON 

INTL  TEL  S  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 


KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  HFARHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.O  DRAWER  oa ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DASIAC 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  S  SPACE  CO.,  INC 
P.O.  SOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 

OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01CY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 

OICY  ATTN  RICHARD  G.  JOHNSON. 

DEPT  52-12 

OICY  ATTN  J.S.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

OICY  ATTN  R.  HEFFNER 


J AYCOR 

11011  TORREYANA  ROAD  M.I.T.  LINCOLN  LABORATORY 

P.O.  BOX  85154  P.O.  BOX  73 

SAN  DIEGO,  CA  92138  LEXINGTON,  HA  02173 

OICY  ATTN  J.L.  SPERLING  OICY  ATTN  DAVID  M.  TOWLE 

OICY  ATTN  L.  LOUGHLIN 

JOHNS  HOPKINS  UNIVERSITY  OICY  ATTN  0.  CLARK 

APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 
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ONMEL  DOUGLAS  CORPORATION 
11  BOLSA  AVENUE 
ITINGTON  BEACH,  CA  92647 
IlCr  ATTN  N.  HARRIS 
nCT  ATTN  J.  MOUUE 
)1CY  ATTN  GEORGE  MR07 
nCY  ATTN  W.  OLSON 
)1CY  ATTN  R.W.  HALPRIN 
)1CY  ATTN  TECHNICAL 

LIBRARY  SERVICES 

SSION  RESEARCH  CORPORATION 

5  STATE  STREET 

VTA  BARBARA,  CA  93101 


ITRE  CORP 

ESTGATE  RESEARCH  PARK 
J20  DOLLY  MADISON  BLVD 
CLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
01CY  ATTN  U.  FOSTER 

ACI FIC-SIERRA  RESEARCH  CORP 
2340  SANTA  MONICA  BLVD. 

OS  ANGELES,  CA  90025 
01CY  ATTN  E-C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  •.OORESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01CY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 


aicY 

ATTN 

P.  FISCHER 

BELLEVUE,  WA 

98009 

Q1CY 

ATTN 

W.F.  CREVIER 

OICY  ATTN 

E.J.  FREMOUW 

01  CY 

ATTN 

STEVEN  L.  GUTSCHE 

01CY 

ATTN 

R.  BOGUSCH 

PHYSICAL  DYNAMICS,  INC. 

01  CY 

ATTN 

R.  HENDRICK 

P.O.  BOX  10367 

OICY 

ATTN 

RALPH  KILB 

OAKLAND,  CA 

94610 

01CY 

ATTN 

DAVE  SOWLE 

ATTN  A.  THOMSON 

OICY 

ATTN 

F.  FAJEN 

OICY 

ATTN 

M.  SCHEIBE 

R  8  D  ASSOCIATES 

01  CY 

ATTN 

CONRAD  L.  LONGMIRE 

P.O.  BOX  9695 

OICY 

ATTN 

B.  WHITE 

MARINA  DEL  REY,  CA  90291 

OICY 

ATTN 

R.  STAGAT 

OICY  ATTN 

FORREST  GILMORE 

OICY  ATTN 

WILLIAM  B.  WRIGHT 

SSION 

RESEARCH  CORP. 

OICY  ATTN 

WILLIAM  J.  KARIAS 

'20  RANDOLPH 

1  ROAD,  S.E. 

OICY  ATTN 

H.  ORY 

.BUQUERQUE, 

NEW  MEXICO  87106 

01 CY  ATTN 

C.  MACDONALD 

OICY 

R.  STELLINGWERF 

OICY  ATTN 

R.  TURCO 

01  CY 

M.  ALME 

OICY  ATTN 

L.  DeRAND 

OICY 

L.  WRIGHT 

OICY  ATJN 

W.  TSAI 

ITRE  CORPORATION,  THE 

RAND  CORPORATION,  THE 

.0.  BOX  208 

1700  MAIN  STREET 

iOFORO 

,  MA 

01730 

SANTA  MONICA, 

CA  90406 

OICY 

ATTN 

JOHN  MORGANSTERN 

OICY  ATTN 

CULLEN  CRAIN 

OICY 

ATTN 

G.  HARDING 

OICY  ATTN 

ED  BEDROZIAN 

01  CY 

ATTN 

C.E.  CALLAHAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01CY  ATTN  VINCE  TRAPANI 


i 


SCIENCE 

APPLI 

CATIONS,  INC. 

VISIDYNE 

1150  PROSPECT 

PLAZA 

SOUTH  BEDFORD  STREET 

LA  JOLLA 

,  CA 

92037 

BURLINGTON,  MASS  01803 

01  CY 

ATTN 

LEWIS  M.  LINSON 

OICY  ATTN  W.  REIDY 

01  CY 

ATTN 

DANIEL  A.  HAMLIN 

01CY  ATTN  J.  CARPENTER 

01  CY 

ATTN 

E.  FRIEMAN 

01CY  ATTN  C.  HUMPHREY 

01  CY 

ATTN 

E.A.  STRAKER 

01  CY 

ATTN 

CURTIS  A.  SMITH 

UNIVERSITY  OF  PITTSBURGH 

PITTSBURGH,  PA  15213 

SCIENCE 

APPLI 

CATIONS,  INC 

OICY  ATTN;  N.  ZABUSKY 

1710  GOOORIDGE  DR. 

MCLEAN, 

VA  22102 

DIRECTOR  OF  RESEARCH 

01  CY 

J.  COCKAYNE 

U.S.  NAVAL  ACADEMY 

01  CY 

E.  HYMAN 

ANNAPOLIS,  MD  21402 

SRI  INTERNATIONAL 

02CY 

333  RAVENSWOOO  AVENUE 

MENLO  PARK,  CA  94025 

01  CY 

ATTN 

J.  CASPER 

01  CY 

ATTN 

DONALD  NEILSON 

01CY 

ATTN 

ALAN  BURNS 

01  CY 

ATTN 

G.  SMITH 

01CY 

ATTN 

R.  TSUNOOA 

01  CY 

ATTN 

DAVID  A.  JOHNSON 

01  CY 

ATTN 

WALTER  G.  CHESNUT 

01  CY 

ATTN 

CHARLES  L.  RINO 

01  CY 

ATTN 

WALTER  JAYE 

01  CY 

ATTN 

J.  VICKREY 

01  CY 

ATTN 

RAY  L.  LEADABRAND 

01CY 

ATTN 

'G.  CARPENTER 

01  CY 

ATTN 

G.  PRICE 

01  CY 

ATTN 

R.  LIVINGSTON 

01  CY 

ATTN 

V.  GONZALES 

01  CY 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 


75  WIGGINS  AVENUE 

BEDFORD, 

MA  01730 

01CY 

ATTN 

W. 

P. 

BOQUIST 

TOYON  RESEARCH 

CO 

. 

P.O.  Box 

6890 

SANTA  BARBARA 

CA 

93111 

OICY 

ATTN 

JOHN 

ISE,  JR. 

OICY 

ATTN 

JO 

EL 

GARBARINO 

TRW  OEFE 

NSE  & 

SPA 

CE  SYS  GROU 

ONE  SPAC 

E  PAR 

K 

REDONDO 

BEACH 

r 

CA 

90278 

01  CY 

ATTN 

R. 

K 

.  PLEBUCH 

01  CY 

ATTN 

S. 

ALTSCHULER 

01  CY 

ATTN 

D. 

DEE 

01  CY 

ATTN 

D/ 

S 

TOCKWELL 

SNTF/ 

1575 
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